NASA 





'lit!*' 



NASA 
Reference 
Publication 
1 233, Vol. 4 

1990 


Satellite Radar 
Altimetry Over Ice 


Volume 4 — Users' Guide for 
Antarctica Elevation Data 
From Seasat 


NASA 

National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Division 


H. Jay Zwally 

Goddard Space Flight Center 

Greenbelt, Maryland 

Judith A. Major and Anita C. Brenner 
ST Systems Corporation 
Lanham, Maryland 

Robert A. Bindschadler 
Goddard Space Flight Center 
Greenbelt, Maryland 

Thomas V. Martin 
Van Martin Consulting, Inc. 

Rockville, Maryland 



TABLE OF CONTENTS 


Pag e 



PREFACE 

V 

1.0 

INTRODUCTION 

1 

2.0 

GEO-REFERENCED DATA BASE 

5 

3.0 

POLAR STEREOGRAPHIC ELEVATION GRID 

27 

4.0 

GEOID GRID 

31 


TABLES 

33 


APPENDIX 

47 


REFERENCES 

54 


iii 

on (I mrrnm/^Pf bwnk 

PRECEDING PAGE BLANK NOT FILMED 





PREFACE 


A gridded surface-elevation data set and a geo-referenced data base for the Seasat radar 
altimeter data over Greenland are described in this volume. It is intended to be a "user’s guide" 
to accompany the data provided to data centers and other users. The grid points are on a polar 
stereogr?phic projection with a nominal spacing of 20 km. The gridded elevations are derived 
from the elevation data in the geo-referenced data base by a weighted fitting of a surface in the 
neighborhood of each grid point. The gridded elevations are useful for the creating of large-scale 
contour maps, and the geo-referenced data base is useful for regridding, creating smaller-scale 
contour maps, and examinating individual elevation measurements in specific geographic areas. 
Tape formats are described, and a FORTRAN program for reading the data tape is listed and 
provided on the tape. For more details of the data processing procedures and corrections that 
were derived and applied to the data, see Volume 1 of this series. 
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SECTION 1.0 
INTRODUCTION 


This volume is the second in a series documenting the data-processing methods and ice 
data products derived from satellite radar altimeter measurements over the ice sheets of 
Greenland and Antarctica and surrounding sea ice. A gridded elevation data set and a 
geo-referenced data base for the Seasat radar altimeter data over Greenland are described in 
this volume. It is intended to be a "user’s guide" to accompany the data provided to data centers 
and distributed to various users on a magnetic tape. The gridded elevations are useful for 
creating large-scale contour maps, and the geo-referenced data base is useful for regridding, 
creating smaller-scale contour maps, and examinating individual elevation measurements in 
specific geographic areas. For more details of the data processing procedures and corrections 
that were derived and applied to the data, see Volume 1 of this series. 

The gridded elevations are on a polar stereographic projection with a nominal spacing 
of 20 km between grid points. The gridded elevation value for each grid point is derived from 
the geo-referenced data base by a weighted fitting of a biquadratic function (or a bilinear function) 
to the elevation data that fall within a certain radius of the grid location. The geo-referenced 
data base contains surface elevations ordered in geographic bins. 

The input Seasat radar altimeter data, in the form of Geophysical Data Records (GDR’s) 
and Sensor Data Records (SDR’s) produced by NASA’s Seasat project at the Jet Propulsion 
Laboratory, was obtained from the NOAA Environmental Satellite Data and Information Service 
(EDIS) archive on about 1000 magnetic tapes. Development of the data processing methods, the 
production of higher-level geophysical data products, and analysis and evaluation of the data have 
been supported at the Goddard Space Flight Center by funding for research and data analysis, 
provided primarily by NASA’s Ocean Processes Program and by the Climate program. Computer 
programming and technical assistance has been provided by the EG&G Washington Analytical 
Services Center, Inc. through December 1988 and by ST Systems Corporation since January 
1989. Numerous other individuals have provided valuable assistance. 

Results have been reported in refeered scientific literature (e.g., Brenner et al., 1983; 
Martin et al., 1983; Zwally et al., 1983; Thomas et al. ( 1983; and Gundestrup et al., 1986). In 
addition, elevation data in various forms have been provided to other scientists and placed in the 
National Snow and Ice Data Center (NSIDC) and the National Space Science Data Center (NSSDC). 
The purpose of this series of reports is to document technical details and provide guidance to 
users of the ice data products. 
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While all reasonable quality-control efforts have been made to eliminate erroneous data, 
some data of questionable quality is likely to have persisted, particularly in the lower-level data 
products. Users should apply normal standards of scientific caution in their use of the data. 

The current list of reports is: 

Satellite Radar Altimetry over Ice.Volume 1: Processing and Corrections of Seasat Data 

over Greenland . July 1989. NASA Ref. Publ. . 

Satellite Radar Altimetry over Ice. Volume 2: User's Guide for Greenland Elevation Data 

from Seasat , July 1989. This volume. 

Satellite Radar Altimetry over Ice. Volume 4: User's Guide for Antarctica Elevation Data 

from Seasat . July 1989. NASA Ref. Publ. 

Volume 3 will be the Antarctic equivalent of Volume 1. Additional volumes will include 
descriptions of the data sets being produced by NASA from the radar altimeter data acquired by 
the U.S. Navy’s GEOSAT, using methods similar to those for the Seasat data. 

The Seasat spacecraft {e.g., Lame and Born, 1982 and Lame et aL, 1980} was launched 
in late June 1978, and during its brief, 110-day lifetime, collected 90 days of nearly continuous 
radar altimeter data from July 9 through October 10 between the latitudes of 72°S and 72°N. 
Although designed only for measurements over water, the Seasat radar altimeter (MacArthur, 
1978; Tapley et al. t 1982; and Townsend, 1980), acquired more than 600,000 useful altimeter 
range measurements over the continental ice sheets of Greenland and Antarctica. 

Over sloping and undulating surfaces, such as ice covered land, or surfaces with highly 
variable reflecting characteristics, such as in regions of sea ice, the range to the surface and the 
characteristics of the received radar pulse changed faster than the response capability of the 
altimeter electronics. Consequently, it has been necessary to correct each range value for lags of 
the altimeter range servo-tracking circuitry by a procedure called retracking (Martin et aL, 1983). 
The retracking correction typically had a mean value of 4- 1.4 m as applied to the surface 
elevation, a standard deviation of 2.9 m, and maximum and minimum values of ± 15 m. In 
addition, the pulse-limited footprint (1.6 km minimum diameter), which was located near the 
satellite nadir point over the relatively flat ocean, was in general located anywhere within the 
beam-limited footprint (22 km in diameter) over sloping surfaces. The resulting slope-induced 
error, which was nearly 80 m over slopes of 0.8 degree, can be partially corrected using the 
procedures described in Brenner et aL, 1983. Corrections are also made for errors in orbit 
determination, atmospheric propagation path-length variations, and earth and ocean tides. 
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Elevation measurements were obtained at 0.1-sec intervals, corresponding to 662 m 
intervals along the subsatellite ground track. The precision of the corrected range measurements 
is about 1.6 m overall with a minimum of about 0.25 m in the smoothest regions of the ice sheets 
(Zwally et al., 1983). The 5- to 10-cm precision over the ocean is for 1-sec data averages. The 
absolute accuracy of the elevations is primarily determined by the limitations on the correction 
methods for the slope-induced errors and by uncertainties in the geoid reference level. 

The principal ice data sets produced and retained are: 

Level 4: Contour maps and gridded elevations with respect to earth ellipsoid and sea 

level (e.g., this volume and Volume 2). 

Level 3: Geo-referenced data base including all individual elevation measurements 

(including time, latitude/longitude positions, and slope-correction estimates) accessible 
by geographic cells (e.g., this volume and Volume 2). 

Level 2: Ice Data Records (IDR’s). Orbital -format data records including altimeter 

parameters, corrected elevations, latitude/longitude positions, AGC, applied corrections, 
retracking beta parameters, and estimates of along-track and cross-track slope corrections. 
(See Volumes 1 and 3.) 

Level 1: Waveform Data Records (WDR’s). Orbital -format data records including waveform 
amplitudes by gate, ranges, AGC, and latitude/longitude positions. (See Volumes 1 and 
3.) 

Altimeter Sensor Data Records (SDR’s) 

The magnetic tape with the gridded elevation and geo-referenced data base was generated 
on an IBM 3081. The data sets are contained on several files of this unlabeled, 6250-bpi tape. 
The geo-referenced data base is written on the first two files (see Tables 2 and 3). A FORTRAN 
program, which can be used to unload and read the data base on files 1 and 2 on the IBM 3081, 
is written on file 3 in ASCII. A listing of this source may also be found in the Appendix. The 
elevation grid over the Antarctic ice sheet is written on files 4 and 5 (see Tables 4 and 5). It is 
important to note that the elevations in the data base are relative to the ellipsoid, while the grid 
elevations are relative to sea level. The Goddard Earth Model 10-B (GEM10-B) geoid grid which 
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was used to obtain elevations relative to sea level is written on files 6 and 7 (see Tables 6 and 
7). File 8 contains detailed information in ASCII concerning the location of various revs and the 
number of points in the gco-rcfcrenccd data base in order of geographical area. File 9 contains 
a narrative description of the tape including the version number, dates of the data, and specific 
information on the sources used to reduce the raw data to surface measurements and grid values 
(see Table 8). All files except files 3, 8 and 9 have been written in IBM binary integer format. 
Blocksizcs vary for each file and arc given in the tables of file descriptions. 
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SECTION 2.0 

GEO-REFERENCED DATA BASE 


The geo-referenced data base contains surface height measurements derived from Scasat 
altimetry data, ordered by geographic areas or "bins". The distribution of the Scasat data used 
in the data base is shown in Figure 1. Figures 2 through 9 show the configuration of the 36,180 
bins in the vicinity of Antarctica. Bin sizes vary in order to compensate for the higher data 
density near Seasafs maximum extent in latitude. Each bin is assigned a number starting with 
1 at the Greenwich meridian. Bin numbers increment first in an eastward direction and then 
from south to north. The starting bin numbers for each row are indicated in the left margin of 
the map in the figures, while the number of data points is printed within the appropriate bin. 
Table 1 is a sample page of the information contained on file 8 of the tape and summarizes the 
number of points and the rev numbers found in each bin, along with the geo- referenced 
coordinates of the southwestern-most corner of the bin. Only bins which contain data are listed 
on file 8 of the tape. The table is written as fixed block in ASCII with a record length of 132 
bytes and blocksize of 19,008 bytes. 

The geo-referenced data base is structured such that the data are ordered first by bin 
number and then by time within each bin. Each data point within each bin contains information 
relating to the position, rev number, surface height relative to the ellipsoid, slope correction and 
orbit adjustment. Corrections which have been applied to the surface elevations are indicated by 
the altimetry data status word in the data base header record found on file 1 of the accompanying 
tape. A detailed explanation of all the corrections may be found in Reference 5. The orbit 
adjustment has been applied to the surface elevation when it was available. Records for which 
the orbit adjustment is unavailable (as indicated by a value of -999999999 in bytes 21-24) have 
the unadjusted surface elevation in bytes 9-12. The user should be aware that using all the 
surface elevation values without checking if the orbit adjustment was valid or not will result in 
an inconsistent data set. The slope correction values are supplied on each data record but have 
not been applied to the surface elevation. The orbit adjustment will improve on the radial 
accuracy of the orbit. The slope correction will compensate for the fact that the original altimeter 
height is measured to the closest point within the radar beam, which is not necessarily the 
subsatellite point. When the slope correction is unavailable a value of -999999999 is placed in 
that field. In order to obtain a slope corrected surface elevation relative to the ellipsoid the 
following algorithm would be used: 


^ h cor - h db ■ ^ h slope 


( 1 ) 
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Figure 5. Seasat Antarctic Geo-referenced Data Base Configuration: 120°E-160°E. 
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Figure 6. Seasat Antarctic Geo-referenced Data Base Configuration: 160°E-200°E. 
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Figure 7. Seasat Antarctic Geo-referenced Data Base Configuration: 240°E-280°E. 
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Figure 9. Seasat Antarctic Geo-referenced Data Base Configuration: 320°E-360°E. 
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where 

^ h cor is surface elevation with the slope correction applied 

h db is the surface elevation in the data base 

AH SL ope * s slop 0 correction. 

To remove the orbit adjustment the following algorithm should be used: 

Ah unadj = H db + ^ h orb (2) 


where 

h unadj is surface elevation without the orbit adjustment 

h db is the surface elevation in the data base 

^ORB is the orbit adjustment 

The data base is designed to be used on a direct-access device, so that data from one or 
several bins may be accessed without the need to read all the records prior to the location 
desired. This is achieved by dividing the data base into three sections. 

The first section of the data base, a header which may be found on file 1 of the 
accompanying tape, gives a summary of its configuration: the locations of the corners of the data 
base, the number of latitude rows, the width in degrees of each of these rows, and the number 
of longitude divisions in each row. These pieces of information give the layout of the data base, 
as depicted in Figures 2 through 9. Information pertaining to the size of the data base, the 
starting record of the bin directory, and the corrections applied to the data are also contained in 
this header. 

Following the header and contained on file 2 of the tape are the altimetry data ordered 
by bin number and within each bin by time. The altimetry data are subdivided into two 
subgroups for each bin which contains data. The first subgroup consists of one logical record 
which indicates the number of data points contained in the bin. The second subgroup consists 
of the actual altimetry data (position, rev number, surface height, orbit adjustment and slope 
correction), with each record corresponding to a data point. 

The final section, a bin directory also contained on file 2, starts at the logical record 
indicated In the data base header. The directory contains an entry for each bin, and starting 
with the first bin, indicates the record number in the data base at which the start of the data 
from a particular bin may be found. Bins which contain no data have a zero entered in the 
directory. Tables 2 and 3 summarize the structure of the data base header and data base in 
greater detail. 
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The data base may be used to locate data within any desired area. The following example 
demonstrates how this may be done. The limits of the desired area are used in conjunction with 
the header Information to determine exactly which bin numbers contain the data. Using the 
southernmost latitude of the desired area along with the width of the latitude rows, establishes 
the southernmost row which contains the data. Longitude limits of the desired area are then 
checked in conjunction with the size and location of the longitude divisions in that row. When 
the longitude limit of the desired area for that latitude group is exceeded, the process starts again 
with the next latitude row to the north. These steps are repeated until the northernmost 
boundary limit of the desired area is reached. 

Equipped with the bin numbers which contain the data, the directory, which gives the 
logical record on the direct-access disk at which each bin begins, is read. If the directory value 
for the bin is non-zero, this logical record is then read to determine the number of records which 
follow and are contained in the same bin. The subsequent data is then read for each bin. 

Software has been developed for use on the IBM 3081 which reads the geo-referenced 
data base on the first two files. A program which reads and prints out the contents of every bin 
given the southeastern and northwestern latitude-longitude limits of a desired area is listed in the 
Appendix and may be found on file 3 of the accompanying tape. The file is in ASCII, is fixed- 
blocked with a record length of 80 bytes, and is blocked at 3,200 bytes. Latitudes should be 
input in degrees North and longitudes in positive degrees East. The subroutine RANDRD along 
with its entry point RANDWR read and write one logical record of data, respectively, utilizing a 
system supplied direct access FORTRAN I/O package which includes DREAD and DWRITE. The 
entry points BLKRD and BLKWR read and write blocks of data at a time. 
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SECTION 3.0 

POLAR STEREOGRAPHIC ELEVATION GRID 


A grid was generated using the corrected and adjusted surface elevations in the 
geo-referenced data base after applying the slope correction. Data for which either the orbit 
adjustment or slope correction were unavailable were not used for the grid. Elevations in the grid 
were obtained by taking data located in the vicinity of each grid point and fitting them to a 
bilinear or biquadratic surface to determine the surface height at the grid point. GEM 10-B geoid 
values were subtracted from the elevations so that they are relative to sea level (see Section 4.0). 

The accompanying grid was generated in a tangent polar stereographic projection where 
the plane of projection is located at the geographic South Pole (the projection latitude) and is 
normal to the earth’s axis. Figure 10a depicts the concept behind this type of projection. A 
straight line is drawn from the North Pole (pole of projection), through a point of the earth’s 
surface, 9, to the projection plane which is tangential to the South Pole. The projection plane 
is in turn divided into square grids from the pole to the Equator with the South Pole at the 
center. Three projection parameters define the size of the plane and the orientation of the plane 
and grid size: 

S - a conversion factor from half-inch grids at the projection latitude to the desired 
grid size; 

<f > p - the minimum latitude extent of the map perimeter for the projection latitude 
located at the North Pole; the maximum latitude extent for the projection latitude 
located at the South Pole; 

G - the Greenwich orientation in degrees. 

In the case of Antarctica, where 20 km grid cells were decided as being optimum for the data 
distribution, values of S=1.65, <f > p = -50°, and G=270° were chosen. 

These three parameters are sufficient to define a grid of the southern hemisphere, from 
the South Pole to -5 <f latitude where the number of grids of desired size from the pole to the 
Equator may be represented by: 

D = — = — . (3) 

S x 10 6 
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where R is the radius of the earth measured in one half-inch grid cells and was chosen to be 
consistent with polar stereographic projections described in other documents. 

The integer number of grids of desired size from the pole to the map perimeter is: 

N — D x tan — ^ 

2 

The grid, defined by I and J axes, with the origin in the upper left corner (see Figure 
10b), represents the coordinates of the South Pole as: 

Ip = N + 1 (5) 

Jp = N + 1 . 

Any point with latitude <f> and longitude A which is located in the southern hemisphere 
south of 4> p is positioned at the following I, J coordinates: 

I - INT[d x A x cos(X) + Ip 4- 0.5] ( 6 ) 

J = INT[d x sin(X) + Jp + 0.5] 

where 

d is D x tan 90 ~ IIeL 
2 

X is A + G 

A is +1 if 4> p > 0 

A is -1 if <f> p < 0 . 

The included grid was generated such that smoothed heights relative to the ellipsoid are 
located at each of the I, J coordinates within Antarctica. Grid locations outside Antarctica and 
any undefined points within Antarctica are indicated by a —100000000. Figure 1 1 depicts a 
topographic map obtained from this grid contoured at 100-m intervals. 

File 4 of the accompanying tape contains a grid header which gives information defining 
the polar stereographic projection used. File 5 contains the grid points which were obtained 
using either a biquadratic or bilinear fit. Details concerning the grldding procedure may be found 
in Reference 5. Data are stored on file 5 such that the information for ten grid points is 
contained in one block of data. The order of grid points is first from decreasing to increasing I, 
then from decreasing to increasing J. Tables 4 and 5 give detailed description of files 4 and 5. 
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Figure 11. Seasat Antarctic Topographic Map Contoured in 100-Meter Intervals. 
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SECTION 4.0 
GEOID GRID 


In order to obtain ice sheet elevations relative to sea level, the geoid was subtracted from 
each grid elevation. Geoid values were bilinearly interpolated from the lxl-degree GEM10-B 
(Goddard Earth Model 10-B) geoid grid. Figure 12 shows a contour of this geoid in the vicinity 
of Antarctica. Files 6 and 7 of the accompanying tape contain the header information and GEM 
10-B lxl-degree grid. Tables 6 and 7 give detailed descriptions of these flies. Any values of the 
geoid not located on the map in Figure 1 1 are set equal to -100000000 in the grid. 
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Sample of ASCII Contents on File 8 Record Format: Fixed Block 

Record Length: 132 Bytes 

Blocksize: 19008 Bytes 



Table 2. Seasat Geo-referenced Data Base Header Description 


GEO-REFERENCED DATA BASE HEADER RECORD 

Record Format: One logical record corresponds to one physical record 



Blocksize: 

424 Bytes 

Bvtes 

FORTRAN 

Variable 

Tvoe 

DescriDtion 

1-4 

1*4 

Number of latitude rows in the data base (49). 

5-8 

1*4 

Northwestern-most latitude of data base in degrees North (x 10 5 ) 
(-6299999). 

9-12 

1*4 

Northwestern-most longitude of data base in degrees East (x 10 ) 
(0). 

13-16 

1*4 

Southeastern-most latitude of data base in degrees North (x 10 5 ) 
(-7209998). 

17-20 

1*4 

Southeastern-most longitude of data base in degrees East (x 10 ) 
(36000000). 

21-216 

1*4 

Width of each latitude row in degrees (x 10 5 ), starting with the 
southernmost row. This is dimensioned by the number of latitude 
rows in the data base. 

217-412 

1*4 

The number of longitude divisions in each latitude row, starting 
with the southernmost row. This is dimensioned by the number 
of latitude rows in the data base. 

413-416 

1*4 

Logical record in data base at which directory starts. 

417-420 

1*4 

Size of the data base, including the directory, in blocks. 

421-424 

1*4 

Status word for altimetry data. 


0 31 
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Table 2 . 
(421-424 cont.) 


Seasat Geo- referenced Data Base Header Description (Cont.) 
Bits Value Description 


0-23 0 

24 1 
0 

25 1 
0 

26 1 
0 

27 1 
0 

28 1 
0 

29 1 
0 

30 1 
0 

31 1 


0 


Unused 

Slope correction applied 
Slope correction not applied 
Orbit adjustment applied 
Orbit adjustment not applied 
Solid tides removed 
Solid tides not removed 
Retracking correction applied 
Retracking correction not applied 
Center of gravity bias applied 
Center of gravity bias not applied 
Tropospheric correction applied 
Tropospheric correction not applied 
Ionospheric correction applied 
Ionospheric correction not applied 
Time bias applied 
Time bias not applied 
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Table 3. Seasat Geo-referenced Data Base Description 


FILE 2: GEO-REFERENCED DATA BASE 

Record Format: 595 logical records correspond to one physical record 



Blocksize: 

19040 Bytes 

Subgroup 1: 

One logical record for each bin containing data 

Bvtes 

FORTRAN 

Variable 

Type 

Description 

1-4 

1*4 

Indicates the number of logical records which follow which are 
located in the bin 

5-32 


Unused 

Subgroup 2: 

One logical record for each data point in the bin 


FORTRAN 

Variable 


Bvtes Type Description 


1-4 

f -T r 

1*4 

North latitude of datum point In degrees (x 10 6 ) 

5-8 

1*4 

East longitude of datum point in degrees (x 10 ) 

9-12 

1*4 

Surface height relative to the ellipsoid in cm. 

13-16 

1*4 

Height sigma, arbitrary value of 1.0 m used (x 10 5 ) 

17-18 

1*2 

Rev number 

19-20 

1*2 

Used for temporary flags when gridding the data 

21-24 

1*4 

Orbit adjustment in meters (x 10 5 ) (-999999999 if unavailable) 

25-28 

1*4 

RMS of orbit adjustment in meters (x 10 5 ) (-999999999 if 
unavailable) 

29-32 

1*4 

Slope correction in meters (x 10 5 ) (-999999999 if unavailable) 

NOTE: 

Subgroups 1 and 2 are 

repeated for as many bins with data. 
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Table 3. Seasat Geo-referenced Data Base Description (Cont.) 


Subgroup 3: Directory 


Bvtes 

FORTRAN 

Variable 

Type 

Description 


1*4 

Record number at which data for bin 1 starts 

5-8 

1*4 

Record number at which data for bin 2 starts 

9-12 

1*4 

Record number at which data for bin 3 starts 

13-16 

1*4 

Record number at which data for bln 4 starts 

17-20 

1*4 

Record number at which data for bin 5 starts 

21-24 

1*4 

Record number at which data for bin 6 starts 

25-28 

1*4 

Record number at which data for bin 7 starts 

29-32 

1*4 

Record number at which data for bin 8 starts 


NOTE: The directory contains as many 32-byte logical records as necessary to designate 

the record locations of all bins. ° 
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Table 4. Elevation Grid Header Description 


FILE 4: ELEVATION GRID HEADER RECORD 

Record Format: One logical record corresponds to one physical record 



Blocksize: 

80 Bytes 

Bvtes 

FORTRAN 

Variable 

Type 

Descriotion 

1-4 

1*4 

Number of latitude increments in the grid for a non-polar 
stereographic grid (294) 

5-8 

1*4 

Number of longitude Increments in the grid for a non-polar 
stereographic grid (294) 

9-12 

1*4 

Starting north latitude of grid in degrees North (x 10 6 ) (this will 
be approximate for a polar stereographic grid) (-73000000) 

13-16 

1*4 

Starting east longitude of grid in degrees East (x 10 6 ) (this will 
be approximate for a polar stereographic grid) (-1800000000) 

17-20 

1*4 

Ending north latitude of grid in degrees North (x 10 6 ) (this will 
be approximate for a polar stereographic grid) (-63000000) 

21-24 

1*4 

Ending east longitude of grid in degrees East (x 10 6 ) (this will be 
approximate for a polar stereographic grid) (1800000000) 

25-28 

1*4 

Status word for data used to generate grid. A zero in any bit 


position indicates that the correction is not applied. 


1 

0 


31 

Bits 

Value 

Descriotion 

0-23 


Unused 

24 

1 

Slope correction applied 


0 

Slope correction not applied 

25 

1 

Orbit adjustment applied 


0 

Orbit adjustment not applied 

26 

1 

Solid tides removed 


0 

Solid tides not removed 

27 

1 

Retracking correction applied 


0 

Retracking correction not applied 

28 

1 

Center of gravity bias applied 


0 

Center of gravity bias not applied 

29 

1 

Tropospheric correction applied 


0 

Tropospheric correction not applied 

30 

1 

Ionospheric correction applied 


0 

Ionospheric correction not applied 

31 

1 

Time bias applied 


0 

Time bias not applied 
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Table 4. Elevation Grid Header Description (Cont.) 


Bvtes 

29-32 

33-36 

37-40 

41-44 

45-48 

49-52 

53-56 

57-60 

61-64 

65-68 

69-72 

73-76 

77-80 


FORTRAN 

Variable 

Type 

1*4 

1*4 

1*4 

1*4 

1*4 

1*4 

1*4 

1*4 

1*4 

1*4 

1*4 

1*4 

1*4 


Description 

Polar stereographic grid size conversion and scaling factor from 
half-inch grids on projection plane to the desired grid size (x 10 6 ) 
(1650000) 

The number of grids of desired size from the pole to the equator 
based on the grid size conversion and scaling factor (x 10 6 ) 
(608754894) 

Latitude of the map perimeter in degrees North (x 10 6 ) 
(-500000000) 

Greenwich orientation in degrees (x 10 6 ) (270000000) 

Polar stereographic switch (1) 

=0, grid has constant Increment in latitude and longitude 
= 1 , grid is in polar stereographic projection 

Number of Taxis divisions to the extent of the map perimeter 
(445) 


Number of J-axis divisions to the extent of the map perimeter 
(445) 

J coordinate of the projected pole (223) 

I coordinate of the projected pole (223) 

Minimum J index of the grid (76) 

Maximum J index of the grid (369) 

Minimum I index of the grid (76) 

Maximum I index of the grid (369) 
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Table 5. Elevation Grid Description 


FILE 5: ELEVATION GRID DATA RECORD 

Record Format: 10 logical records correspond to one physical record 



Blocksize: 

1800 Bytes 

Bvtes 

FORTRAN 

Variable 

Tvoe 

Descriotion 

1-4 

1*4 

Condition number of the matrix used in the least squares solution 
to the function (x 10 6 ) 

5-8 

1*4 

Capsize in degrees latitude - radius from grid location defining 
area from which data was used to define grid (x 10 ) 

9-12 

1*4 

North latitude of grid point in degrees (x 10 ) 

13-16 

1*4 

g 

East longitude of grid point in degrees (x 10 ) 

17-20 

1*4 

Height values of the grid at location relative to sea level in meters 
(x 10 5 ) 

21-24 

1*4 

Number of data values that were used to calculate grid value 

25-28 

1*4 

Number of parameters used to define function, NPT, (equals 0, 
3, or 6) 

29-52 

1*4 

Six gridding function coefficients. If NPT is < 6 then the rest of 
the coefficients are initialized to zero, (x 10 5 ) 

53-76 

1*4 

Set of null coefficients associated with any negligible singular 
values (see SVD reference). If NPT Is < 6 then rest of coefficients 
are initialized to zero (x 10 6 ) 

77-80 

1*4 

g 

Distance in km from grid locations to closest data point (x 10 ) 

81-84 

1*4 

North latitude of closest data point to grid location in degrees 
(x 10 6 ) 

85-88 

1*4 

East longitude of closest data point to grid location in degrees 
(x 10 6 ) 

89-92 

1*4 

Height associated with closest data point to grid location in 
meters (x 10 5 ) 

93-96 

1*4 

Standard deviation of the data with respect to the gridding 
function in meters (x 10 6 ) 

97-180 

1*4 

Correlation matrix from solution. This is a symmetrical 6X6 
matrix so only the upper triangular portion is stored. The order 
of storage is elements 1-6 are the first row elements, 7-11 
columns 2-6 of second row etc. (x 10 ) 

NOTE: 

Ten of the above-mentioned 180-byte logical records make up one block of data. 
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FILE 6: 


Bvtes 

1-4 

5-8 

9-12 

13-16 

17-20 

21-24 

25-28 

29-32 

33-36 

37-40 

41-44 

45-48 

49-52 


Table 6. Geold Grid Header Description 
GEOID GRID HEADER RECORD 

Record Format: One logical record corresponds to one physical record 

Blocksize: 80 Bytes 


FORTRAN 

Variable 

Type Description 

IM Number of latitude Increments in the grid for a non-polar 

stereographic grid (16) 

** 4 Number of longitude increments in the grid for a non-polar 

stereographic grid (361) 

** 4 Starting north latitude of grid in degrees (x 10 6 ) (this will be 

approximate for a polar stereographic grid) (-750000000) 

I>4 Starting east longitude of grid in degrees (x 10 6 ) (this will be 

approximate for a polar stereographic grid) (0) 

I<4 Ending north latitude of grid in degrees North (x 10 6 ) (this will 

be approximate for a polar stereographic grid) (-60000000) 

r * 4 Ending east longitude of grid in degrees East (x 10 6 ) (this will be 

approximate for a polar stereographic grid) (3600000000) 

1*4 Unused 


l* 4 Polar stereographic grid size conversion and scaling factor from 

half-inch grids on projection plane to the desired grid size (x 10 6 ) 
(1650000) 6 

** 4 The number of grids of desired size from the pole to the equator 

based on the grid size conversion and scaling factor (x 10 6 ) 
(608754894) 

I* 4 Latitude of the map perimeter in degrees North (x 10 6 ) 

(-500000000) 

I<4 Greenwich orientation in degrees East (x 10 6 ) (2700000000) 

!* 4 Polar stereographic switch (1) 

=0, grid has constant increment in latitude and longitude 
= 1, grid is in polar stereographic projection 

I>4 Number of I-axis divisions to the extent of the map perimeter 

(445) 
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Table 6. Geoid Grid Header Description (Cont.) 


FORTRAN 

Variable 


Bvtes 

Type..- 

DescriDtion 

53-56 

1*4 

Number of J-axis divisions to the extent of the map perimeter 
(445) 

57-60 

1*4 

J coordinate of the projected pole (223) 

61-64 

1*4 

I coordinate of the projected pole (223) 

65-68 

1*4 

Minimum J index of the grid (59) 

69-72 

1*4 

Maximum J index of the grid (386) 

73-76 

1*4 

Minimum I index of the grid (59) 

77-80 

1*4 

Maximum I index of the grid (386) 
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Table 7. Geoid Grid Description 


FILE 7: GEOID GRID RECORD 

Record Format: 200 logical records correspond to one physical record 



Blocksize: 

2400 Bytes 

Bvtes 

FORTRAN 

Variable 

Tvoe 

DescriDtion 

1-4 

1*4 

North latitude of grid point in degrees North (x 10 6 ) 

5-8 

1*4 

East longitude of grid point in degrees East (x 10 6 ) 

9-12 

1*4 

Value of geoid in meters (x 10 5 ) 


NOTE: Two hundred of the above-mentioned 12-byte logical records make up one block of 
data. 
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Table 8. Narrative Description of Tape 


FILE 9: NARRATIVE DESCRIPTION 

Record Format: One logical record corresponds to one physical record 

Blocksize: 80 Bytes 


Record 1: 


Bvtes 

FORTRAN 

Variable 

Tvdc 

Description 

1-8 

A*8 

Satellite name 

9-11 

A*3 

Version number of this tape 

12-17 

A*6 

Date of release of this data (YYMMDD) 

18-80 

A*63 

Differences between this data release and previous versions 


Record 2: 

FORTRAN 

Variable 

Bytes Type 

1-80 A* 1 


Description 

Continuation of differences between this data release and previous 
versions 


Record 3: 


Bvtes 

FORTRAN 

Variable 

Type 

Description 

1-10 

AMO 

Orbit used to compute surface elevations 

11-20 

AMO 

Gcoid used to compute elevations relative to sea level 

21-40 

A*20 

Source of ocean tides 

41-60 

A*20 

Source of troposphere correction 

61-80 

A*20 

Source of ionosphere correction 
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Table 8. Narrative Description of Tape (Cont.) 


Record 4: 


FORTRAN 

Variable 


Bvtes 

Tvoe 

DescriDtion 

1-10 

AMO 

Surface used for orbit adjustment 

11-40 

A*30 

NASA publication number and title for documentation describing 



data on tape 

41-46 

A*6 

Beginning day of data (YYMMDD) 

47-52 

A*6 

End day of data (YYMMDD) 

53-80 

A*28 

Name of mission 
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APPENDIX 
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START 

COL 


1 

7 


Program to Load Data Rase onto Direct Access Device and Read Data for Specific Area 




— 1 — 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

7 

7 

7 

7 

6 

7 

7 

7 

6 

7 

6 

7 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


-+ 7 +- 


FUNCTION: GIVEN THE LATITUDE AND LONGITUDE LIMITS OF A RECTANGULAR ar fa 
THIS PROGRAM DETERMINES WHICH DATA BASE BINS ARE CONTATNFn tkj* 
THE AREA AND PRINTS OUT THE CONTENTS OF THOSE BINS N 


INPUT - UNIT 5 


UNIT 12 
UNIT 20 
UNIT 21 


ooI^I“5 0UTHEASTERN no RTH LATITUDE OF DESIRED AREA 

EAST LONGITUDE OF DESIRED AREA 
GNWLAT-NORTHWESTERN NORTH LATITUDE OF DESIRFD ARFA 

SCRATCH N DIsl( W FnR E nsT« 4 |Lt 0 ?P ITUDE 0F DESIRED "REA 
Mt"S! 5 “hea!>ER <15? FILE " 3350 ” ISK) 

DATA BASE ENTIRES AND DIRECTORY (2ND FILE OF TAPE) 


□ T0 READ antarctic D.B. PGMR - J.A. 

REVISED 12/17/86 FOR WDCG TAPE PGMR - J.A 


MAJOR 

MAJOR 


C 

C 


DOUBLE PRECISION XLAT, XLON 
INTEGERX2 IPASS,IFLAGS 
LOGICALX1 LKICK 

DIMENSION ZNWION(2),ZSELON<2) 0) 

DIMENSION JJBIN(IOO), LKICK (100) 

COMMON/RANBLK/LAT,LON,IHT, ISIG, IPASS, IFLAGS, IORB, IORBR, ISLCR, 

COMMON/CHEAD/IG, IWNLAT, IWNLON, ISELAT , ISELON, IWDLAT, LONDIV, IDRREC, 

n?VPxk^f E o« IDRD(1),l - AT) ' fIHEAD(l),IG) 

. DAT N100/100/ 0/,IFIL2/21/ ' I0UT6/6/ ' NBINS/36180/ 'IN5/5/, 

DATA LKICK/100*. FALSE./ 

READ LIMITS OF DESIRED AREA 

5pT?pJTm',^ 0 ?^ n ^ ELAT ' GSEL0N ' GNNLA TAGNWLON 

WRITE(I0UT6, 10003) GSELAT, GSELON, GNWLAT, GNWLON 

READ HEADER FROM DATA BASE TAPE 
READ(IFIL1, 10001) IHEAD 

C0NV WNLAr= P IWNLATSl R E-5 LUES T ° FL0ATING P0INT 

E-5 
E-5 


50 


NNLON = IWNLONX1 

SELAT = ISELATX1 

SELON = ISEL0NX1.E-5 
DO 50 1=1, IG 
WIDLAT ( I ) = IWDLAT(I)*1 .E-5 
CONTINUE 

WRITE(I0UT6, 10000) WNL AT, WNL ON, SELAT, SELON, IG, WIDLAT, LONDIV, 
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Program to Load Data Base onto Direct Access Device and Read Data for Specific Area (Cont.) 


START 

COL + 1 


--+ 7 + 8 


6 

1 

7 

7 

1 

7 

1 

7 

7 

1 

1 

7 

7 

7 

7 

7 

1 

7 

7 

1 

1 

1 

7 

7 

7 

7 

6 

7 

7 

7 

7 

7 

7 

1 

7 

1 

7 

7 

1 

1 

7 

7 

7 

7 

7 

7 

1 

1 

7 

7 

7 

7 

1 

1 


c co«pute 1> the I size C 6f P the G lo»oitude divisions for each bin group 

SIZE(I) I=1 tSELON-HNLON) / LONDIV(I) 

100 CONTINUE oT , c 

WRITE(I0UT6, 10004) SIZE 
C DUMP REMAINING DATA TO DIRECT ACCESS DISK 
REWIND IFIL1 
CALL BLKRDC IFIL2, NBUF) 


200 


READ DIRECTORY AND PRINT OUT ITS CONTENTS 
DO 200 1 = 1 > NBINS 
IREC = (I-l)/8 + IDRREC 
CALL RANDRD 
IWICH = M0D( 1-1 » 8 ) + 1 
IDIRCI) = IDRD(IWICH) 

CONTINUE 

WRITE(IOUT6,2210) 

WRITE(I0UT6,2220) IDIR 


C DETERMINE BIN NUMBER LOCATED AT SOUTHWEST CORNER OF DESIRED AREA 
C AND STORE THE START RECORD OF THAT BIN 

XNWLON = DM0D( GNWL ON+36 0 . DO , 36 0 . DO ) 

fF<ON!ilAT D TsllAV OM 5R‘ 0 G§ELAT 0 GE 0 iNlAT .OR. GNHLON .GE. SELON 
.OR. GSELON . LE . WNLON) GO TO 900 
DIF = GSELAT - SELAT 
ADD = 0.0 
DO 320 IB=1,IG 
ADD = ADD + WIDLAT (IB) 

IBGRP1 = IB „ 

IF(DIF .LT. ADD) GO TO 325 
CONTINUE 
I BGRP 1 = I G 

CLAT1 = ADD - WIDL AT ( IBGRP1 ) 

^Jt^G^WNLAT. CNHUAT = HNLAT 

C CHECK IF GREENWICH SPLITS THE CAPSIZE REGION 

IFCGNWLON . LT T 0 . DO^AND . G GSELON^GE . 0 . DO ) GO TO 327 
IA = 1 

XNWLON 
XSELON 


DATA 112 
DATA 113 


320 

325 


ZNWLON(l) 
ZSELON(l) 
GO TO 328 
REGION SPLIT. 


327 


IA = 2 
ZNWLONC 1 ) 
ZSELON(l) 
ZNWL0NC2) 
ZSEL0NC2) 


.DIVIDE INTO TWO PIECES 

XNWLON 
360. DO 
0 . DO 
XSELON 


C 

c 


LOOP OVER THE AREA(S) CONTAINED BY THE CAPSIZE REGION 


DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 


129 

130 

131 

132 

133 
139 

135 

136 

137 

138 

139 

190 

191 

192 

193 
199 
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Program to Load Data Base onto Direct Access Device and Read Data for Specific Area (Cont ) 

START 

COL — - +— 2 + 3 -—+ , + 5 


328 CONTINUE 

DO 700 M=1 , IA 
CLAT = CLAT1 
IBGRP = IBGRP1 
CSELON = ZSELONCM) 

- GT - selon) cselon=selon 

CNWLON = ZNWLONCM) 

IFCZNWLONCM) .LT. WNLON) CNWLON=WNLON 
DIF = CNWLON - WNLON 
330 IBIN = 0 

IF( IBGRP .EQ. 1) GO TO 345 
IFIN = IBGRP - 1 
DO 340 J = 1 , IFIN 
340 IBIN = IBIN + LONDIVC J) 

345 IBIN = IBIN + DIF/SIZECIBGRP) +1.0 

IF( IBIN .GT. NBINS) GO TO 360 
NBIN = NBIN + 1 
IFCNBIN .GT. N100) GO TO 980 
IIBIN(NBIN) =IBIN 
NRECNOC NBIN )=IDIRC IBIN) 

C LOCA JloS L = d1^S^E?iSgRpI ITUDE GR0UP WHICH ARE C0NTAINED IN DESIRED 
£L0N = WNLON + NLON*SIZECIBGRP) 

350 CLON = CLON + SIZE(IBGRP) 

IFCCLON .GT. CSELON) GO TO 360 

IFCCLON. EQ. SELON) GO TO 360 

IBIN = IBIN + 1 

IF( IBIN .GT. NBINS) GO TO 360 

NBIN = NBIN + 1 

I F( NBIN .GT. N100) GO TO 980 

IIBIN(NBIN) =IBIN 

NRECNOCNBIN)=IDIR(IBIN) 

GO TO 350 

mS R 0 C !ut T ° ClS IF 0UISIDE DESIRED ARE ‘ 

TOTLAT = SELAT + CLAT + .0001 
IBGRP = IBGRP + 1 
IFC IBGRP .GT. IG) GO TO 700 
IFCTOTLAT .LT. CNWLAT) GO TO 330 
700 CONTINUE 

C RE "°g S40 L |=?Tnb?S CURANCES ° F BIN NUMBERS DUE T0 GREENWICH CROSSING 
DO 420 K = 1,'nBIN 
IF( LKICKC I ) ) GO TO 420 
IFC I .EQ. K) GO TO 420 

=r?\t!!uV IBIN<K>) 00 T0 420 

420 CONTINUE 
440 CONTINUE 
NJBIN = 0 
DO 450 1=1, NBIN 
IFC LKICKC I) ) GO TO 450 
NJBIN = NJBIN+1 
JJBINCNJBIN) = IIBINCI) 


DATA 145 
DATA 146 
DATA 147 
DATA 148 
DATA 149 
DATA 150 
DATA 151 
DATA 152 
DATA 153 
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START 
COL - 


— + 1 + 2 + 


3 + 4 + 5 + 6 + 7 


+ - 


1 

7 

1 

1 

1 

7 

7 

7 

7 

7 

7 

1 

1 

7 

7 

7 

7 

1 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

6 

1 

1 

7 

1 

1 

1 

7 

1 

1 

7 

1 

2 

2 

6 

2 

2 

2 

6 

2 

6 

1 

6 

6 

6 

6 


450 ^rxTEC I0UT6 ,30000) NJBIN,CJJBINCJ),J=1,NJBIN) 


C STORE RECORD NUMBER OF BIN TO BE READ; IF ZERO, THEN BIN CONTAINS 
C NO DATA. . .PROCEED TO NEXT BIN 
DO 550 1=1, NBIN 
IFC LKICKC I ) ) GO TO 550 
IREC = NRECNOC I ) 

IFCIREC .NE. 0) GO TO 460 
WRITECI0UT6, 40000) IIBIN(I) 

C FIRST°RECORD°OF BIN TELLS HOW MANY RECORDS IN THE BIN 
460 CONTINUE 

CALL RANDRD 

WRITEC I0UT6, 2400) IIBINCI ) , NRECNOC I) , LAT 
WRITE(IOUT6,2100) 

C PRINT OUT CONTENTS OF ONE ENTIRE BIN 
DO 500 J = 1 , 1 END 
IREC = IREC + 1 
CALL RANDRD 
ORBIT = IORBXl.E-5 
ORBITR = IORBRXl . E-5 
SLPCOR = ISLCRX1.E-5 
SIGH = ISIGXl.E-5 
XL AT = LATXl.D-6 

WRITEC I 0UT6^ 2200) IREC, XL AT , XLON, IHT , I PASS, ORBIT, ORBITR, SLPCOR, 

. SIGH 

500 CONTINUE 
550 CONTINUE 
STOP 
c 

C AREA REQUESTED OUTSIDE OF DATA BASE AREA 
900 WRITECI0UT6, 90000) 

STOP 

C TOO MANY BINS FOR ARRAY WHICH STORES DIRECTORY 
980 WRITECIOUT6, 81000) N100 
STOP 
c 

2200 FORMATCIIO, 2F10 . A, 110, 18, 2F10 . 3, F10 .3, Fll .3) 

2100 F0RMATC/4X, ' IREC ', 5X, ' LAT 1 , 9X, ' LONG' , 7X, ' HEIGHT ', 3X, REV # , 
4X,'0RB ADJ',2X,'0RB RMS' ,2X, 'SLOPE CORR' » AX, 'SIGMA' ) 

2210 FORMATC///' DIRECTORY FOR DATA BASE (IDIR)»'/) 

2300 FORMAT (/// ^DIRECTORY STARTS AT LOGICAL RECORD', 110, 

2AOOFORMAT(/' M FOR BIN^NUMBER' , I6^2X^ ' LOGICAL RECORD STARTS AT', 110, 
2X, ' t RECS IN BIN:', 110) 

10000 FORMAT(//’ HEADER INFORMATION FROM DATA BASE TAPE'/ 

2X,'NW LATITUDE & LONGITUDE (WNLAT, WNLON) i , 2F10 . 9/ 

' 2X *SE LATITDUE 8 LONGITUDE CSELAT , SELON) : • , 2F10 . A/ 

2X 'LAT WIDTH OF EACH OF THE', 14,' BIN GROUPS: ’/ 

6(1X,8F9 .2/), F10.2/2X, 'LONGITUDE DIVISIONS IN EACH BIN GROUP', 
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START 

COL + 


2 + 3 + 4 + 5 + 6 + 7 


6 

6 

6 

1 

1 

1 

6 

1 

6 

1 

6 

1 

1 

6 

1 

6 

7 

1 

1 

1 

7 

1 

1 

1 

1 

1 

1 

1 

1 

1 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

4 

7 

4 

7 

4 

4 

7 

1 

7 

7 

7 

7 


10001 

10002 

10003 

10004 

30000 

40000 

81000 

90000 


‘ , ^ONDIV) : V6(1X,8I6/),I6/2X, f NUMBER OF BUFFERS : M6 , 3X, 

1 LOGICAL RECORD AT WHICH DIRECTORY STARTS: ' , 110// 

1 ALTIMETRY DATA STATUS W0RD',Z1Q//) 

FORMAT ( 120A4) 

FORMAT ( 4F10 . 3 ) 

FORMAT (// * SOUTHEASTERN CORNER OF DESIRED AREA ( LAT, LON) : 1 , 2F10 3 
• nD ^! t^ 0 ? THWESTERN C0RNER 0F D^IRED AREA ( LAT, LON) , 2F10 ! 3/) 
F0 R 7(1X / 8F9 I 3/)/) EACH BIN GR0UP AS DETERMINED FROM HEADER:*/ 

FORMATC//^ 1/ jqqq j G qz ) G M6 ,* BINS CONTAIN DATA IN THE DESIRED* 

FORMAT (/ * BIN *,I6,» CONTAINS NO DATA*) 

FORMAT (// ' *x NUMBER OF BINS CONTAINED IN AREA EXCEEDS * , 15, 

INCREASE SIZE OF NRECNO ARRAY'/ 1 PROGRAM TERMINATING') 
FORMAT (/// * AREA SELECTED IS OUTSIDE DATA BASE AREA'/ 

' PROGRAM TERMINATING') 

END 


CC 

CC 

CC 


C 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE RANDRD 


FUNCTION; 


THIS ROUTINE USES THE DIRECT ACCESS I/O PACKAGE CDAIO) 

AND ENTRY POINTS TO ACHIEVE SEVERAL THINGS: 

ll damSmS ~ THE IREC TH LOGICAL RECORD FROM DISK 

£f^™ R ” WRITES THE IREC TH LOGICAL RECORD TO DISK 

3) BLKRD - TRANSFERS THE DESIRED PHYSICAL RECORDS 

OR BLOCKS FROM UNIT NIN TO DISK 

4) BLKWR - TRANSFERS THE DESIRED PHYSICAL RECORDS 
OR BLOCKS FROM DISK TO UNIT NOUT 


C0MM0N/FERMSG/IMESC26) 

REAL A,BUF(8,595),TBUF(8,595) 
INTEGER DISK 


C0MM0N/RANBLK/AC8), IREC 
EQUIVALENCE (TBUF( 1,1), BUF(1, 1 ) ) 

DATA NREC/595/, NW0RDS/8/, KBUF/ 1 /, DISK/12/ 

LOGICAL WSWTCH 

K=1 


GO TO 10 
ENTRY RANDWR 
WSWTCH= . TRUE . 
K = 2 


10 I BUF= ( IREC-1 )/NREC+l 

IF CIBUF.NE.KBUF) GO TO (50,90),K 
20 JREC=IREC-CIBUF-1)XNREC 
GO TO (30,70),K 
30 DO 40 I =1 , NWORDS 
40 A(I)=BUF(I, JREC) 

RETURN 

50 IF (WSWTCH) CALL DWRITEC DISK, KBUF, BUF) 
KBUF=IBUF 
WSWTCH=. FALSE. 


CALL DREADCDISK, KBUF, BUF, 860) 
GO TO 20 


+ 8 
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START 

COL 


Program to Load Data Base onto Direct Access Device and Read Data for Specific Area (Cont.) 


+ 1 + 2 + 3 + A + 5 + 6- + - 7 

60 PRINT 1000, IMES 
STOP 

70 DO 80 1=1 , NW0RDS 

80 BUF( I , JREC)=A( I ) 

RETURN 

90 CALL DWRITEC DISK, KBUF, BUF) 

KBUF=IBUF 
GO TO 20 

1000 FORMAT (1H0,Z8,I6,20A4,4(1X,Z8)) 

ENTRY BLKRDC NIN, NBUF) 

KBUF=NBUF 
DO 205 J=1 , NBUF 
READ(NIN,1200,END=204) TBUF 
205 CALL DWRITEC DISK, J , BUF) 

20A WRITEC 6 ,1100) J 

1100 FORMAT ( 1 NUMBER OF BUFFERS READ ’,15) 

1200 FORMAT (18(25 5A4),170A4) 

RETURN 

ENTRY BLKWRC NOUT , NBUF) 

DO 200 J = 1 , NBUF 

CALL DREADCDISK, J, BUF, 860) 

200 WRITEC NOUT ) BUF 
RETURN 
END 
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